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We study the temperature-dependent dephasing rate of excitons in chains of chromophores, ac- 
counting for scattering on static disorder as well as acoustic phonons in the host matrix. From this 
we find a powerlaw temperature dependence of the absorption line width, in excellent quantitative 
agreement with experiments on dye aggregates. We also propose a relation between the line width 
and the exciton coherence length imposed by the phonons. The results indicate that the much de- 
bated steep rise of the fluorescence lifetime of pseudo-isocyanine aggregates above 40 K results from 
the fact that this coherence length drops below the localization length imposed by static disorder. 
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A wide variety exists of materials consisting of in- 
teracting chromophores; examples are molecular aggre- 
gates, conjugated polymers, natural light-harvesting sys- 
tems, and arrays of quantum dots. Lately, such multi- 
chromophore systems have received much attention for 
their collective optical properties; 1 ! 2 ! 3 ! 4 efficient and con- 
trolled energy transport processes^ the occurrence of 
strong quantum entanglement of collective chromophore 
states§*£ and their possible application in schemes for 
quantum computation^^ All these properties sensitively 
depend on the temporal and spatial energy and phase 
relaxation of the chromophores' collective excitations, 
which make such relaxation processes a topic of broad 
interest and importance. The nature of the collective ex- 
citations and their dynamics are governed by the compli- 
cated interplay of their scattering on static disorder (e.g., 
inhomogeneity in the excitation energy of the individual 
chromophores) and on dynamic degrees of freedom (vi- 
brations). A variety of experimental techniques are used 
to investigate the resulting intricate dynamics. The the- 
oretical study of the combined effects of static disorder 
and a heat bath on exciton decoherence and optical re- 
sponse remains relatively unexplored. 9 

Recently, a model of weakly-localized Frenkel excitons 
coupled to acoustic phonons in the host medium proved 
to yield an excellent description of the nonmonotonous 
temperature dependence of the fluorescence Stokes shift 
as well as the linear temperature scaling of the radiative 
lifetime of aggregates of the dye 3,3'-bis(sulfopropyl)-5,5'- 
dichloro-9-ethylthiacarbo-cyanine (THIATS) measured 
between and 100 Ki4 More direct insight into the 
interplay of static disorder and homogeneous dephas- 
ing, and thus a more stringent test for the model, is 
provided by temperature dependent absorption, photon 
echo, and hole burning measurements. Such measure- 
ments, over the interval 0-300 K, have been performed 
for the prototypical J-aggregates of the dye pseudo- 
isocyanine (piC)ii2iiAii2iiiii4 No unique interpretation of 
the thermal dephasing and line broadening in multichro- 
mophore systems has emerged from these studies: acti- 
vated processes involving several optical phonons of the 



aggregate have been suggested as mechanism; 11 ! 12 ! 13 but 
also power-law broadening due to host vibrations^ An- 
other unsolved issue is the temperature dependence of the 
fluorescence lifetime at elevated temperatures. For aggre- 
gates of PIC-Br, its steep rise above 40 K has been only 
explained using a two-dimensional aggregate model, 15 
which contradicts the generally accepted chain-like ge- 
ometry of these systems in solution^ 

In this Letter, we revisit the model of excitons in a lin- 
ear chain with static energy disorder, coupled to acoustic 
phonons in the host. We show that the dephasing rates 
calculated within this model give an excellent explana- 
tion of the temperature dependent absorption measure- 
ments for PIC aggregates. Furthermore, for the strong 
scattering regime we introduce a relation between the 
exciton coherence length and the dephasing rates and 
demonstrate that this explains the steep rise of the fluo- 
rescence lifetime. 

Our model consists of a linear Frenkel chain of N two- 
level chromophores (n — 1,...,N), with uncorrelated 
Gaussian disorder of standard deviation a in the chro- 
mophore energies e n and resonant dipole-dipole interac- 
tions between chromophores n and m given by J nrn = 
— J/\n — m| 3 , J > being the nearest-neighbor inter- 
action. All transition dipoles are parallel. The exci- 
ton eigenstates (labeled v) follow from diagonalizing the 
N x N matrix with the e„ as diagonal elements and the 
J nm as off-diagonal ones. They have energy E v , and site- 
amplitudes ip vn ] the disorder leads to their localization 
on segments of the chain. The model includes on-site 
scattering of excitons on a phonon bath. If the exciton- 
phonon interaction is not too strong, it induces scattering 
from exciton state v to state fi with a rate that may be 
obtained from Fermi's Golden rule (see, e.g. Ref. 0): 

N 

= t(M)gm J2 vW vn ■ (i) 

Here, w M „ =E li - E v and T(uj) = 2nJ2 q \Vq\ 2 S(u - u q ), 
the one-phonon spectral density {uj q is the energy of a 
phonon in mode q and V q characterizes the coupling of 
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this mode to the excitons; we set K = 1). Furthermore, 
G(u) = n(u) if id > and G(ui) = 1 + n(—u) if ui < 
0, with n(uj) = [exp(w/fc B T) - l] -1 , the mean thermal 
occupation number of a phonon mode of energy ui. 

While in Ref . U the were only used to study the in- 
traband exciton relaxation and the resulting fluorescence 
kinetics, they also govern the temperature dependent de- 
phasing of the excitons, and hence their homogeneous 
line widths. Specifically, the thermal dephasing rate of 
the ^th exciton state is given by T v = | and de- 

pends on temperature through the n(u>). In contrast to 
the fluorescence kinetics, the temperature dependence of 
the dephasing rates and absorption line width sensitively 
depends on the form assumed for the spectral density. 
Here, we will assume that J-(to) = Wq(ui / J) 3 . This form 
is appropriate for acoustic phonons, for which the density 
of states (DOS) scales as ui 2 (Debye behavior) and \V q \ 2 
scales like lo; the remaining factor Wq is a free parameter 
in the model. As we will see, the lo 3 scaling gives a natu- 
ral explanation for the line width measurements on linear 
dye aggregates. Similar good agreement between theory 
and experiment is found if we allow for fluctuations of 
J- (to) around an average u> 3 scaling, but the agreement is 
lost by changing to a power different from 3i£ 

Using the above model, we simulated the absorption 
spectrum for an ensemble of chains. The expression for 
this spectrum reads: 



A(E) 




r„ + 7„/2 



7T {E-E„) 2 + {T„+ lv /2) 2 
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where F v — (YlnVvn) 2 is the dimensionless oscillator 
strength of the vth exciton state, "f v = joF v is its radia- 
tive rate (70 is the emission rate of a single chromophorc) 
and the angular brackets denote averaging over the static 
disorder in the e n . The resulting line width A(T) is de- 
termined as the full width at half maximum (FWHM) 
of the simulated spectrum. We used chain lengths of 
250-1000 chromophores and averaged over 1000-400,000 
disorder realizations, depending on parameter choices. 

In the absence of disorder and for long chains with 
nearest- neighbor interactions only, we found that A(T) oc 
T 3 ' 5 (neglecting the radiative contribution, which is al- 
lowed unless T is very low). The explanation is simple: 
for cr = 0, the lowest (y = 1) exciton state dominates the 
spectrum. The above scaling then simply results from 
summing the W^^=\ over all higher lying states /1. The 
lo 3 dependence of T{oj), combined with the lj~ x / 2 de- 
pendence of the exciton DOS near the bottom of the 
one-dimensional band, upon integration yields the T 3 ' 5 
behavior. It turns out that the long-range dipolar inter- 
actions slightly change the power from 3.5 to 3.85ii£ 

If we include disorder, the spectral width at zero 
temperature, A(0), is dominated by the inhomogeneous 
width, which for linear aggregates scales like cr 4 / 3 ^ From 
extensive simulations we have found that in the presence 
of disorder the thermal broadening follows a power-law: 
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FIG. 1: (a) Calculated temperature-dependent width of 
the absorption spectrum (symbols) and corresponding fits 
to Eq. (lines) for six values of the degree of disorder 
a/ J = 0.05, 0.1, 0.2, 0.3, 0.4, 0.5 (bottom to top). The fit pa- 
rameters (a,p) take the values (1.1, 4.0), (1.2, 4.2), (1.3, 4.3), 
(1.2, 4.3), (1.2, 4.3), and (1.0, 4.1), respectively. W /J = 25 
and 'yo/J = 1.5 x 10 . (b) Comparison of theory (curves) 
and experiment (symbols) for the absorption line width as a 
function of temperature for aggregates of PIC-C1 (triangles) 
and PIC-F (squares). The model parameters are given in 
the text. The inset compares the line shape calculated for 
PIC-C1 aggregates (curves) to experiment (symbols) for three 
temperatures. We note that the thermal shift of the absorp- 
tion band is not accounted for in our theory; the calculated 
spectra were shifted by hand to facilitate comparison of the 
line shape. Experimental data were taken from Ref. Il4 



where a and p depend only weakly on the parameters a 
and Wq. The scaling relation Eq. © implies that, al- 
though the spectrum results from a distribution of exci- 
ton states with different dephasing rates, the total width 
may be interpreted as the sum of an inhomogeneous 
width, A(0), and a dynamic contribution. For six a/ J 
values in the interval [0.05, 0.5] (taking Wo/ J — 25), our 
simulated data are presented as symbols in Fig.^a), to- 
gether with fits to Eq. • The fit parameters are given in 
the caption. The power p is seen to lie in a narrow range, 
p = 4.2±0.2, which is close to the disorder-free value 3.85. 
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The reason for this small effect of disorder is that the dy- 
namic width is governed by scattering from the optically 
dominant DOS tail states to high-lying intra-band states, 
whose DOS is modified only slightly by disorder. The 
small increase of p primarily results from the possibility 
to scatter further downward in the tail of the disordered 
exciton band. 

A scaling relation like Eq. JSJ was also found experi- 
mentally by Renge and Wild 14 for aggregates of PIC with 
counter ions Cl~ and F - . Fits of their data, presented 
in Fig. n^b), demonstrate good agreement between the- 
ory and experiment. These fits were obtained using the 
accepted values of J = 600 cm -1 and 70 = 2.7 x 10 8 
s _1 (= 1.5 x 10~ 5 J) for PIC and taking the fit param- 
eters a/ J = 0.128 and W /J = 16.0 for PIC-C1 and 
a I J = 0.249 and W /J = 15.9 for PIC-F. In the inset 
we show (for PIC-C1) that not only the width, but also 
the experimental line shape as a function of temperature 
is recovered very well by our model. 

We next address PIC-Br aggregates, for which the ab- 
sorption line width as well as the fluorescence lifetime 
were measured up to 190 K (triangles and squares, re- 
spectively, in Fig. We first discuss the width, for 
which it turned out that our model also for this counter- 
ion gives a very good fit to experiment (solid line in Fig. [2 
fit parameters <r/J = 0.135, Wo/J = 32.5). It should 
be noted that the same data have previously been fitted 
assuming activated behavior, supposedly resulting from 
scattering on optical phonons of the aggregated In or- 
der to describe the strong increase of the line width for 
large temperatures, such a fit requires both a high ac- 
tivation energy (several hundreds of cm -1 ) and a large 
pre-exponential factor. The pre-factor is directly related 
to the exciton- vibration coupling constant. A simple per- 
turbative treatment shows that the thus estimated cou- 
pling constant should be so large (few thousand cm _1 's) 
that a strong polaron effect is to be expected. This is 
inconsistent with the generally accepted excitonic nature 
of the optical excitations in cyanine aggregates. 
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FIG. 2: Calculated and measured line width and fluorescence 
lifetime for aggregates of PIC-Br as a function of temperature. 
Line width: triangles (measured—) and solid (calculated); 
lifetime: squares (measured^) and dashed (calculated). 



At low temperature, the fluorescence lifetime is deter- 
mined by a competition between spontaneous emission 
and intraband relaxation of the weakly localized exci- 
tons, and may be analyzed using a master equation for 
these populations. 4 While successful for THIATS aggre- 
gates up to 100 K, this approach turns out not to explain 
the steep rise of the measured lifetime for PIC-Br above 
about 40 K (cf. Fig. [SJ. In fact, using the above fit pa- 
rameters, the model only yields a doubling of the lifetime 
from to 150 K, which is an order of magnitude too small. 
Other approaches have seen similar problems in explain- 
ing the measurements jiSiisi the closest (far from perfect) 
agreement has been found using a two-dimensional ag- 
gregate modeliiii Here we argue that the problem arises 
from phonon-induced decoherence of the excitons inside 
their localization segments. 

If we consider weak exciton-phonon scattering in ho- 
mogeneous chains (a = 0), the exciton coherence length 
(or mean- free path), A™ h , is much larger than its wave- 
length. In other words, the uncertainty in the exciton's 
quasi- momentum, 5K ~ , is much smaller than 

the quasi-momentum itself. Due to the quadratic disper- 
sion relation near the band bottom, = 2 J — JK 2 , 
8K results in an energy uncertainty 6E^ — 2JK5K ~ 
JK/N^. Identifying 5E^ with the exciton's homoge- 
neous width r M , we thus obtain an estimate for the co- 
herence length: A^ oh ~ JK/T^. In the strong scattering 
regime, A™ h becomes smaller than the exciton wave- 
length and the quasi-momentum is not a good quantum 
number anymore: SK ~ K. Hence SE^ ~ J(SK) 2 ~ 
J/(N™ h ) 2 , leading to the estimate N™ h - (J/rj 1 / 2 . 

We now address this issue for the localized band edge 
states in the presence of disorder. At low temperature 
[ksT <C A(Q)]i the coherence length A™ h associated 
with the exciton-phonon scattering is much larger than 
the localization length A^ 00 ; hence the latter plays the 
role of the coherence length and determines the opti- 
cal dynamics. The main consequence of the exciton- 
phonon scattering in this regime is the intraband relax- 
ation (equilibration) described above already^ Upon in- 
creasing T, A™ h decreases and at some temperature T 
becomes comparable with N l ° c . In other words, excitons 
start to scatter within localization segments. For T > To, 
the coherence length A^ oh , rather than the localization 
length, governs the optical dynamics. It is to be stressed 
that when T surpasses To, the optically dominant states 
immediately fall in the strong scattering regime, because 
they resemble standing waves without nodes with a typ- 
ical "wavelength" on the order of A^ oc i2£ 

The cross-over temperature To separates two different 
regimes for the fluorescence lifetime as a function of tem- 
perature. This lifetime is inversely proportional to the 
number of coherently emitting chromophores. At low 
temperature, this number should be identified with the 
typical localization length of the optically dominant ex- 
citon states, N*; at high temperatures, it is to be iden- 
tified with the T-dependent coherence length, A coh (T). 
Denoting To = I/70 as the single-chromophore lifetime, 
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the essential behavior of the lifetime is thus given by: 

f r /N* , T < T , 

r(T) = { (4) 

[ r /iV coh (T) = Toby^iTyj , T > T , 

where in the estimate for the coherence length we used 
for the typical dephasing rate the effective homogeneous 
width 2T(T) = A(T) - A(0) and b is a parameter of 
order unity. Here, the value of To is determined from the 
condition by/T(T )/J = l/N*. 

We used this approach to analyze the fluorescence data 
for PIC-Br. We took N* = 54, as follows from the low- 
temperature lifetime r(0) = 68 ps in combination with 
To = 3.7 nsjiSii^, and which also is consistent with the 
participation ratioi^ of 58 calculated for a/J = 0.135. 
The resulting very good fit, using the A(T) calculated 
above, is seen as dashed line in Fig. [5] Here we used 
b = 2.2 as only fit parameter, which gave To = 44 K. 

Two comments are in place. First, in the low- 
temperature regime our fit does not account for pop- 
ulation relaxation, which for growing T increases the 
lifetime from the value tq/N* due to thermal popu- 
lation of dark states. Using the master-equation ap- 
proach of Ref. |4j, we found that this effect for PIC-Br 
leads to a small (~ 20 ps) linear growth of the life- 
time between T = and T = T). Second, a scaling 
of the lifetime with the line width was also suggested 
by Feldmann et al~ for quantum wells. Translating 
their arguments to a one-dimensional system, one obtains 
t(T) = (tq /tt) \J A(T) / J, which in the high-temperature 
limit agrees with our result, except that the pre- factor 
b is reduced considerably. The reason for this difference 



is that the argument of Ref. based on counting the 
number of states within the absorption width, does not 
relate to the coherence length, but rather to population 
redistribution. 

In conclusion, a model that considers weakly-localized 
excitons which scatter on acoustic host phonons, yields 
power-law thermal line broadening that is in excellent 
agreement with experiments on dye aggregates. In addi- 
tion, a new method to relate the exciton coherence length 
within the strong-scattering regime to the homogeneous 
line width, gives good agreement with temperature de- 
pendent fluorescence lifetime data. Our work shows that, 
in contrast to earlier claims, a one-dimensional model 
for PIC aggregates suffices to explain all experimental 
data. We believe that both the on-site (diagonal) dis- 
order and the one-phonon scattering on host vibrations 
with a Debye-like spectral density (oc w 3 ) are the essential 
and physically sound ingredients to reach this conclusion. 
It is known that off-diagonal disorder yields a symmetric 
low-temperature absorption profile^ while two-phonon 
scattering (including pure dephasing) results in power- 
law broadening with a power that is about twice as large 
as the values of p obtained for one-phonon scattering^ 
Both facts are inconsistent with experimental data. Fi- 
nally, off-site scattering on acoustic phonons is also char- 
acterized by a spectral density proportional to uj 3 (see, 
e.g. the work cited in Ref. Il9|) : therefore its inclusion 
would not change the power p — 4.2, but only alter the 
microscopic meaning of the phenomenological scattering 
parameter Wo to a sum of two (on-site and off-site) scat- 
tering parameters. 
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